Numerical solution of perturbed Kepler problem using a split operator technique

An efficient geometric integrator is proposed for solving the perturbed Kepler motion. This method is stable and accurate over long integration time, which makes it appropriate for treating problems in astrophysics, like solar system simulations, and atomic and molecular physics, like classical simulations of highly excited atoms in external fields. The key idea is to decompose the hamiltonian in solvable parts and propagate the system according to each term. Two case studies, the Kepler atom in an uniform field and in a monochromatic field, are presented and the errors are analyzed.Comment: 17 pages, 5 figures, submitted to the Journal of Computational Physic

Quantum dynamics of a high-finesse optical cavity coupled with a thin semi-transparent membrane

We study the quantum dynamics of the cavity optomechanical system formed by a Fabry-Perot cavity with a thin vibrating membrane at its center. We first derive the general multimode Hamiltonian describing the radiation pressure interaction between the cavity modes and the vibrational modes of the membrane. We then restrict the analysis to the standard case of a single cavity mode interacting with a single mechanical resonator and we determine to what extent optical absorption by the membrane hinder reaching a quantum regime for the cavity-membrane system. We show that membrane absorption does not pose serious limitations and that one can simultaneously achieve ground state cooling of a vibrational mode of the membrane and stationary optomechanical entanglement with state-of-the-art apparatuses.Comment: 14 pages, 7 figure

Quantum dynamics of a vibrational mode of a membrane within an optical cavity

Optomechanical systems are a promising candidate for the implementation of quantum interfaces for storing and redistributing quantum information. Here we focus on the case of a high-finesse optical cavity with a thin vibrating semitransparent membrane in the middle. We show that robust and stationary optomechanical entanglement could be achieved in the system, even in the presence of nonnegligible optical absorption in the membrane. We also present some preliminary experimental data showing radiation-pressure induced optical bistability.Comment: 6 pages, 2 figures. Work presented at the conference QCMC 2010 held on 19-23 July 2010 at the University of Queensland, Brisbane, Australi

TALON - The Telescope Alert Operation Network System: Intelligent Linking of Distributed Autonomous Robotic Telescopes

The internet has brought about great change in the astronomical community, but this interconnectivity is just starting to be exploited for use in instrumentation. Utilizing the internet for communicating between distributed astronomical systems is still in its infancy, but it already shows great potential. Here we present an example of a distributed network of telescopes that performs more efficiently in synchronous operation than as individual instruments. RAPid Telescopes for Optical Response (RAPTOR) is a system of telescopes at LANL that has intelligent intercommunication, combined with wide-field optics, temporal monitoring software, and deep-field follow-up capability all working in closed-loop real-time operation. The Telescope ALert Operations Network (TALON) is a network server that allows intercommunication of alert triggers from external and internal resources and controls the distribution of these to each of the telescopes on the network. TALON is designed to grow, allowing any number of telescopes to be linked together and communicate. Coupled with an intelligent alert client at each telescope, it can analyze and respond to each distributed TALON alert based on the telescopes needs and schedule.Comment: Presentation at SPIE 2004, Glasgow, Scotland (UK

SkyDOT (Sky Database for Objects in the Time Domain): A Virtual Observatory for Variability Studies at LANL

The mining of Virtual Observatories (VOs) is becoming a powerful new method for discovery in astronomy. Here we report on the development of SkyDOT (Sky Database for Objects in the Time domain), a new Virtual Observatory, which is dedicated to the study of sky variability. The site will confederate a number of massive variability surveys and enable exploration of the time domain in astronomy. We discuss the architecture of the database and the functionality of the user interface. An important aspect of SkyDOT is that it is continuously updated in near real time so that users can access new observations in a timely manner. The site will also utilize high level machine learning tools that will allow sophisticated mining of the archive. Another key feature is the real time data stream provided by RAPTOR (RAPid Telescopes for Optical Response), a new sky monitoring experiment under construction at Los Alamos National Laboratory (LANL).Comment: to appear in SPIE proceedings vol. 4846, 11 pages, 5 figure

Efficient numerical diagonalization of hermitian 3x3 matrices

A very common problem in science is the numerical diagonalization of symmetric or hermitian 3x3 matrices. Since standard "black box" packages may be too inefficient if the number of matrices is large, we study several alternatives. We consider optimized implementations of the Jacobi, QL, and Cuppen algorithms and compare them with an analytical method relying on Cardano's formula for the eigenvalues and on vector cross products for the eigenvectors. Jacobi is the most accurate, but also the slowest method, while QL and Cuppen are good general purpose algorithms. The analytical algorithm outperforms the others by more than a factor of 2, but becomes inaccurate or may even fail completely if the matrix entries differ greatly in magnitude. This can mostly be circumvented by using a hybrid method, which falls back to QL if conditions are such that the analytical calculation might become too inaccurate. For all algorithms, we give an overview of the underlying mathematical ideas, and present detailed benchmark results. C and Fortran implementations of our code are available for download from http://www.mpi-hd.mpg.de/~globes/3x3/ .Comment: 13 pages, no figures, new hybrid algorithm added, matches published version, typo in Eq. (39) corrected; software library available at http://www.mpi-hd.mpg.de/~globes/3x3

Optomechanical sideband cooling of a thin membrane within a cavity

We present an experimental study of dynamical back-action cooling of the fundamental vibrational mode of a thin semitransparent membrane placed within a high-finesse optical cavity. We study how the radiation pressure interaction modifies the mechanical response of the vibrational mode, and the experimental results are in agreement with a Langevin equation description of the coupled dynamics. The experiments are carried out in the resolved sideband regime, and we have observed cooling by a factor 350 We have also observed the mechanical frequency shift associated with the quadratic term in the expansion of the cavity mode frequency versus the effective membrane position, which is typically negligible in other cavity optomechanical devices.Comment: 15 pages, 7 figure

CFD Code Validation against Stratified Air-Water Flow Experimental Data

Pressurized thermal shock (PTS) modelling has been identified as one of the most important industrial needs related to nuclear reactor safety. A severe PTS scenario limiting the reactor pressure vessel (RPV) lifetime is the cold water emergency core cooling (ECC) injection into the cold leg during a loss of coolant accident (LOCA). Since it represents a big challenge for numerical simulations, this scenario was selected within the European Platform for Nuclear Reactor Simulations (NURESIM) Integrated Project as a reference two-phase problem for computational fluid dynamics (CFDs) code validation. This paper presents a CFD analysis of a stratified air-water flow experimental investigation performed at the Institut de Mécanique des Fluides de Toulouse in 1985, which shares some common physical features with the ECC injection in PWR cold leg. Numerical simulations have been carried out with two commercial codes (Fluent and Ansys CFX), and a research code (NEPTUNE CFD). The aim of this work, carried out at the University of Pisa within the NURESIM IP, is to validate the free surface flow model implemented in the codes against experimental data, and to perform code-to-code benchmarking. Obtained results suggest the relevance of three-dimensional effects and stress the importance of a suitable interface drag modelling

DNA-Protein Binding Rates: Bending Fluctuation and Hydrodynamic Coupling Effects

We investigate diffusion-limited reactions between a diffusing particle and a target site on a semiflexible polymer, a key factor determining the kinetics of DNA-protein binding and polymerization of cytoskeletal filaments. Our theory focuses on two competing effects: polymer shape fluctuations, which speed up association, and the hydrodynamic coupling between the diffusing particle and the chain, which slows down association. Polymer bending fluctuations are described using a mean field dynamical theory, while the hydrodynamic coupling between polymer and particle is incorporated through a simple heuristic approximation. Both of these we validate through comparison with Brownian dynamics simulations. Neither of the effects has been fully considered before in the biophysical context, and we show they are necessary to form accurate estimates of reaction processes. The association rate depends on the stiffness of the polymer and the particle size, exhibiting a maximum for intermediate persistence length and a minimum for intermediate particle radius. In the parameter range relevant to DNA-protein binding, the rate increase is up to 100% compared to the Smoluchowski result for simple center-of-mass motion. The quantitative predictions made by the theory can be tested experimentally.Comment: 21 pages, 11 figures, 1 tabl